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Abstract 

The modified factorization technique of a quantum system characterized by position-dependent 
mass Hamiltonian is presented. It has been shown that the singular superpotential defined in 
terms of a mass function and a excited state wave function of a given position-dependent mass 
Hamiltonian can be used to construct non-singular isospectral Hamiltonians. The method has been 
illustrated with the help of a few examples. 
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I. INTRODUCTION 



The one dimensional time independent position dependent mass Schrodinger equation 
(PDMSE) Q 

p — r^P ] +v{x) 



HiP{x) 



^ = E^, (1.1) 



2mo \ m(x) 

[p = ~ih^ being the momentum operator and m(x), the dimensionless form of mass func- 
tion M(x) = mom(x)], has attracted a lot of attention during the past few years due to its 
applications in describing the dynamics of electrons in many condensed-matter systems 
The concept of position dependent mass comes from the effective mass approximation [3| 
which is an useful tool for studying the motion of carrier electrons in pure crystals and also 
for the virtual-crystal approximation in the treatment of homogeneous alloys as well as in 
graded mixed semiconductors. This field has also arisen due to the impressive development 
in crystallographic growth techniques [^] which allow the production of non uniform semi- 
conductor specimen with abrupt heterojunctions (see e.g., [sl] for a review). In different areas 
of possible applications of low dimensional structures as already mentioned, there is need to 
have Hamiltonians with predetermined energy spectrum. For example, in the quantum well 
profile optimization [gI], isospectral potentials are generated through supersymmetric quan- 
tum mechanics. This is necessary because intersubband optical transitions in a quantum 
well may be grossly enhanced by achieving the resonance conditions e.g. appropriate spac- 
ings between the most relevant states and also by tailoring the wave functions so that the 
matrix elements relevant for this particular effect are maximized . Consequently, there has 
been a growing interest in obtaining exactly solvable PDM Hamiltonians through the usage 
of supersymmetric quantum mechanics (SUSYQM) [S-ll], point canonical transformation 
l3-14|. Lie algebra [l5[, Heun equation [l6[ and intertwining method 12, (and refer- 
ences cited therein). However, despite being studied widely, the number of exactly solvable 
PDM Hamiltonian is still rather few. 

Fairly recently, the standard factorization technique [l9l-21| has been modified [22| so that 



the excited state wave function of a given nonsingular constant mass Hamiltonian can be 
used to construct nonsingular isospectral partner Hamiltonian. In general the superpotential 
defined by excited state wave function gives rise to the singular partner Hamiltonians with 



negative energy eigenvalues and usual degeneracy of the energy spectrum is lost [23|, |2J] . In 
ref.0, the nonuniqueness of the factorization of a constant mass Hamiltonian is used to 
make it possible to construct acceptable (nonsingular) isospectral Hamiltonians. This algo- 
rithm is independent of the choice of wavefunctions used to define the superpotential. Since 
exactly solvable PDM Hamiltonians are few in number, it should be interesting therefore, to 
extend the modified factorization technique to handle cases with position-dependent mass. 
It may be noted that because of the presence of a non-constant mass the algorithm in the 
case of PDMSE is different from the constant mass case and here our objective is to exam- 
ine whether or not factorization based on excited states can be used to search for physically 
acceptable Hamiltonian within the context of PDMSE. 

In order to make this report self contained, in section HA we have briefly recalled the no- 
tion of usual factorization technique of a quantum system with position dependent mass. In 
section IIB, we have generalized the usual factorization technique so that excited state wave 
functions can be used to construct new exactly solvable PDM Hamiltonians. Furthermore, 
as an useful application of the proposal, in Section HI, we have considered two examples. 
Section IV is kept for summary. Throughout this report we shall consider h = 2mo = 1. 
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II. FACTORIZATION IN QUANTUM SYSTEM WITH POSITION DEPENDENT 
MASS 



A. Usual factorization 

The minimal version of PDM SUSYQM is based on ground state wave function and the 
charge operators are represented hj Q = A^a- and Q''' = A^a^, where a± = ai ± 102 are 
the combinations of pauh matrices a\^2 and the associated two first order hnear operators 
are given by jsf 

^ , 1 ( 1 



where 'prime' denotes the differentiation with respect to x. The components of the super- 
symmetric Hamiltonian H = diag{H~,H~^) can be factorized in terms of Aq and Aq to 
obtain ^ ^ 

= AUo = + —T- + 2.2 

m ax^ ax 

where without loss of generality we have considered zero ground state energy of Hq and the 
two partner potentials Vq" can be written in terms of superpotential Wq and mass function, 
as 



V,^{x) = W,'t 



Worn' 1 / 1 



mi y/m \ y/m 



(2.3) 



m J V 2 

Clearly, V^^ = ^0 +^^^:^ (^) ' "^^^ ground state wave function ip^ of the Hamiltonian 
Hq determines the superpotential 

Woix) = , (2.4) 

which implies that for a physically acceptable mass function^ the superpotential and hence 
two partner potentials Vq" are free from singularity (provided the given potential is 
without any singularity). It is well established that when SUSY is unbroken then the 
spectrum of the two PDM partner Hamiltonians are degenerate except for the lowest 
energy of Hq. Moreover, the energy eigenvalues and normalized eigenfunctions of the SUSY 
partner Hamiltonians Hq are related by [sl 

^(0)+ _ ^(0)- ^(0)- _ Q 

= (i^Sr)-^ (2-5) 
= (Ef +)-^ n = 0,l,2... 

Summing up, for a well behaved mass function the usual factorization technique for PDMSE 
allows us to construct nonsingular isospectral potentials of a given nonsingular potential 
using equations (12.41) and (12. 3p . It should be noted here that the factorization technique 
mentioned above is not unique and can be extended to use excited state wave functions. 
This generalization is mathematically straightforward but physically very nontrivial since it 



^ The mass function should be positive and without any singularity through out its domain of definition. 
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yields singular partner potentials. In the following we shall show that the excited state of a 
given position dependent mass Hamiltonian can be used to construct nonsingular isospectral 
potentials. 



B. Excited state factorization and its modification 

In order to generalize the usual factorization technique to excited states we define the su- 
perpotential as 



where ip~ is the n-th excited state wave function of an exactly solvable nonsingular Hamilto- 
nian Hq . It is clear that for a well behaved mass function the superpotential Wn is singular 
function but for a suitably chosen singular mass function this becomes nonsingular. Since 
singular mass functions are physically absurd so we can not consider this type of mass func- 
tion. Now subtracting the energy En'' of the excited state from the Hamiltonian so that 
the resulting Hamiltonian can be factorized, we have the form of the two partner potentials 

'2WL 1 f I \" ^^-^^ 



v;r x) = v- x) + 



The singular superpotential Wn contributes nothing new to except a constant energy 
shift to Vq~ . In this case energy eigenvalues of the Hamiltonians H~ and are related 
by eI'^^ = E^^^ — En '' ,k = 0, 1, 2...(fc denotes the energy level and n refers to the nth 
eigenfunction of H~ used in factorization). But the partner potential becomes singular at 
the node(s) of the wave function ip'. This singularity is responsible for both the destruction 
of degeneracy of the spectrum and creation of negative energy state(s). The singularity 
breaks the real axis into more than one disjoint intervals and imposes additional boundary 
conditions on the wave functions. In this case, the supersymmetry operators map square 
integrable functions to states outside the Hilbert space. As a result, some or all of the wave 
functions of the SUSY partner Hamiltonians may not belong to the same Hilbert space of 
square integrable functions and usual proof of degeneracy between the excited states of the 



partner Hamiltonians does not hold [23|, |24 



To overcome this, we require the modification of this excited state factorization technique 
which can be achieved by deforming the operators^ as 



where Wn is given by equation (12. 6p and the unknown function /„(x) will be evaluated 
shortly. Using this operators we first obtain 

H: = AUn = --f-, +—,^ + ynHx) (2.9) 

m dx^ dx 



are obtained from equation (|2.ip after replacing Wq by Wn 
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where 

= v;t + ^+ f 2iyn + ^]fn + fl (2.10) 



At this point we assume 



+ f2w^ + ) + = /3, /3 e M (2.11) 



y/i^ V 2m2 

so that if^ = if.^ + /9 is solvable. This assumption also ensures the non uniqueness of the 
factorization viz. 

A+A- + ^ = A+A-. (2.12) 

For (3 = 0, the Riccati equation (12.111) reduces to the Bernoulli equation whose solution is 
given by 

fn{x) = y A e M. (2.13) 

For nonzero /3 the Riccati equation (12.111) is not always solvable but it can be transformed 
into the following second order linear differential equation 



1 ,,// m' 



m m 
with the help of following transformation 



ml rn2 \/m \ ^/m 



^+ = (2.14) 



fn{x) = —={logXn)\ Xn = ^ (2.15) 

wm Jm 



Using the relation (12. 7p it is clear that the equation (I2.14p is (nearly) isospectral to the 
following solvable equation 

- -r" + —^r' + y,{x)r = (^f - - (2.16) 

m m^ 

where the solutions ip^lx) corresponding to two equations (I2.14p and (I2.16P are related by 

(2.17) 

By changing the order of operations between A^ and A~ and using relations (12. 7p . (12. lip 
we obtain the partner Hamiltonian H~ of as 

7 7, 1 d"^ m' d ~ , ^ , , 

H- = A-A^ = — + — — + V-ix) 2.18 

m ctx^ dx 

where 

V- = V--^ + P = Vi -^k- + /3. (2.19) 



m \ m 



At this point few comments are worth mentioning. For a given nonsingular potential Vq 
and a mass function the new potential V^(= + /?) is singular but the potential V~ has 
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no singularity if fn{x) is nonsingular in which case the corresponding spectrum of the two 
isospectral Hamiltonians H~ and H~ are related by 

E^"^" = 4"^" + /3, A; = 0,1, 2.... (2.20) 

For /3 = 0, the singularity of /„ given in (12.131) can be controlled by the arbitrary constant A 
and in all other cases it depends on the nature of solution of the equation (I2.16p . Moreover, 
the eigenf unctions of the two nonsingular Hamiltonians H~ and H~ are related by 

rk - A-Airk ^ ■ ^ 

except for the zero energy state of H~ . The proof of this results is as follows: we start with 

(2.22) 



— Al^A^{Aj^A:^)iPi^ +J3A„A^^i^ 
= A~A^{H-ij-) + pA'^Atru = (4"^" + muA^^k- 
In the second step we have used the non-uniqueness of the factorization of and in the 

_ (n)— — — 

last step we have used the PDMSE H~ipf, = E^. . Hence, if ipj^ is a eigenfunction of H~ 
corresponding to eigenvalue -^f."^" then A'A^tp^ is a eigenfunction of H~ with eigenvalue 
E^''~ + (3. The second relation of equation (I2.2ip can be proved in a analogous way. The 
wave function ip^i (say) corresponding to zero energy of the Hamiltonian H~ can be obtained 
by directly solving A^ip^, = 0. 

Summing up, excited state wave function of a given position dependent mass Hamiltonian 
can be used to construct nonsingular partner Hamiltonian H~ using equation (I2.19P provided 
the function /„ is nonsingular. This has been made possible in two step factorization. In the 
first step, an intermediate singular Hamiltonian has been created and in the second step 
non-uniqueness of the factorization technique has been used to remove the singularities. 
The operator — y4~74+ is similar to the second order reducible intertwining operator C = 
^('^)c^ ^(^) Kientioned in ref.jl^] which intertwines two Hamiltonians H~ and 

Hn^n^n^AnAnHn, (2.23) 

if one identifies —rj[x) = + ^ and — 7(x) = V~ + fnWn- In this regard the present 



algorithm is equivalent to the confluent second order SUSY transformation [25[ in constant 
mass case but the relation between the two isospectral partner potentials is different in the 
present approach due to the position dependence of the mass function. Moreover, the energy 
levels of the two nonsingular isospectral partner Hamiltonians H~ and H~ are degenerate. 
In figure 1, we have shown schematically the combined action of the operators A~A'^ and 
A~A'^ on the wave function of two isospectral partner Hamiltonians and H~ respectively. 
The operators A^ or ^4+ destroy a node but A~ or A~ create the same in the eigenfunctions 
so that the overall number of nodes remains the same. 



III. ILLUSTRATION 

It is worth mentioning here that the proposed technique mentioned in the preceding 
section is more general and can be applied to any solvable PDM Hamiltonian to obtain 
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FIG. 1: Energy levels of two partner Hamiltonians H~ , H~ and the intermediate Hamiltonian 
. Combined actions of the operator A~A^ and A~A^ are shown. 



some acceptable isospectral Hamiltonians which might be useful in various fields of condensed 
matter physics. For illustration purpose we are considering here two exactly solvable position 
dependent mass Hamiltonians. In example 1 we have considered /3 = while non zero /3 has 
been considered in example 2. 

Example 1 (/3 = 0): We consider the following mass function and potential possessing 



harmonic oscillator like spectra [26 



"^(^) = 1 , \ o (3.1) 



1 + a^x^ 



The bound state solution are given by 



1 — ; 



-1 ^ 2 

fih (ax) 



1 fiy e '\ " J ^ f smh-\ax) \ 

^^=V^W TTT^W^ 'I — ^ — ) 

Ef^~ = 2k + l, A; = 0,1, 2... (3.4) 

where Hk{x) denotes the Hermite polynomial. Now we consider here first excited state 

ip-^ = — {i+x^y/^ ' ~ 0, and a = 1 for which the partner potential V-^ and Vi 

are obtained using equation (12.71) . (I2.13P and (12.191) as 

Vr = (sinh-V))'-^(^)-3 (3.5) 



= _l±^+(sinh-ix)2-3+2v/(r + ^-^ 



4 (sinh ^ x) ' 



0fg(sinh '^)2(2A + Erf(sinh"^x)) -2sinh"^x 

(3.6) 

respectively. The energy eigenvalues of the Hamiltonian are given 

0, 1, 2... The normalized first excited state wave function corresponding to zero energy of iff 
is evaluated by solving Afip^ = as 



V8A2-2 e- 



sinh ^ X 



;i + x2)i [2^/^Xe(^'''^ ' - 2 sinh"^ X + v^e(^'>^h '^)' Erf (sinh"^ x)] ' 



(3.7) 



Hence the new potential v}^^ is nonsingular for A > 1/2. _The lowest energy wavefunction 
corresponding to negative energy —2 is obtained through A'^A'^tp^ as 



^0 



TT 



4[2A + Erf(sinh"^x)]e- 



;i + x2)4 [2v/7rAe(^i°'^ ' - 2 sinh"^ x + v^vr Erf (sinh"^ x)e('^i°t^ '^■)'] 



(3.8) 



In figure 2(a), we have plotted the potential Vf{x), its isospectral partner V{' for two 
parameter values A = 1, .7 and the mass function. In figure 2(b), lowest and first excited 
state wave functions of the Hamiltonian iff for A = 1 have been plotted. 
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FIG. 2: (a) Plots of the given potentials V^' (solid line) and Vf for A = 1 (dashed line), A = .7 
(dotted line) and the mass function m{x). (b) The lowest (solid line) and first excited state (dashed 
line) wave functions of the Hamiltonian iff for A = 1. 



Example 2 (/3 7^ 0): Here we consider the following mass function and the potential 

m(x) = ^sech' (|) (3.9) 

Vo{x) = i^±i_^l!^e" + ^i^^e-"", c>^,a + b-c+^>0, xe (-00, 00), (3.10) 

for which one linearly independent solution of the equation fl2.16p is known to be (see 
appendix of ref.[l3] and [27j ) 

^ (x) = Ci e"^/^(l + e^) 2 ^Fi{ , , c, e M , (3.11) 



8 



where 



P2 = 

p(o)- 



(a + 6)2 - 2c(a + 6 - c + 1) + 4(Ef " - 



(3.12) 



The bound state wave functions of the Hamiltonian Hq are given [17], in terms of Jacobi 
polynomial vj^' ( 



cz/2 

-p(c-l,a+b-c) 

(1 + e^)('^+''+i)/2 " 



1 + 



n 



0,1,2... 



(3.13) 



Now we consider first excited state i.e. n = 1 for which the function can be obtained 
using equations f l2.15p and f l2.17p as 



fn 



where Af 



m dx 



' (,og[4rr(.)])+^^ 



2c 



[a + h-c+ l)e^' 



a + 6 + 3 
1 + 



6 — 1^ +sinh ( 
(3.14) 



m(x),^/' and ^/'i are given by equations fl3.9p . fl3.1ip and 



fl3.13p respectively. In order to obtain new nonsingular potential V^' we have to choose 
a, 6, c, /3 such that has no singularity. The analytical expression of /„ is too involved so 
it is very difficult to find the range of parameters values for which /„ is nonsingular. In 
figure 3(a), we have plotted the potential Vi{x) = Vq~ — Ef'^ and its nonsingular partner 

= V{ —'ifn/ ~ -^1°'' + P particular values of a = 1,_6 = 5, c = 4 and /3 = 1. For 
the same set of parameter values, isospectral partner potential V2 has been drawn in figure 
3(b) by considering second excited state factorization i.e. n = 2. The energy eigenvalues of 
the new Hamiltonians and H2 are given 

respectively. In both cases the wave functions can be obtained using equations (I2.2ip and 

dmi). 
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FIG. 3: Plots of (a) the potential V-[' = Vq — 13 (solid line), its nonsingular partner V-[' (dashed 
line) and the mass function m(x), (b) the potential = Vq — 22 (solid line) and its partner 
(dashed line) for a = 1, 6 = 5, c = 4, /? = 1. 
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IV. SUMMARY 



To conclude, we have generalized the modified factorization technique of a Schrodinger 
Hamiltonian suggested by Berger and Ussembayev j22| to a quantum system characterized 
by position-dependent mass Hamiltonian. This generalization is done in such a way that 
a given mass function and the excited state wave function of a PDM Hamiltonian can be 
used to generate new solvable physical Hamiltonians. The method discussed here is applied 
to a number of exactly solvable PDM Hamiltonians which leads to nontrivial isospectral 
Hamiltonians. In the constant mass limit of the present algorithm, all the recently published 



results [22| on modified factorization can be recovered. 

It would be interesting to study the modified factorization technique for a quantum system 
characterized by broken supersymmetry. 
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